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AN  INTEGRAL  EQUATION  OCCURRING  IN  A  PROBLEM  OF  SUBSIDENCE 

ABSTRACT 

The  linear  boundary  value  problem  that  was  solved  by  numerical 
methods  in  BRL  Report  No.  1131  has  been  reformulated  as  a  linear 
integral  equation  and  solved  by  Laplace  transform  methods. 
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INTRODUCTION 


As  a  result  of  mining  operations  cavities  are  produced  into  which 
the  overburden  may  eventually  collapse.  Following  collapse  there  will 
be  a  period  during  which  the  fallen  material  subsides .  It  would  be 


extremely  difficult  to  treat  this  motion  in  a  way  that  takes  into 
account  all  the  geometrical  and  physical  peculiarities  of  a  particular 
situation.  However,  Leser  and  Jenike  ^  have  proposed  a  semi-erpirical 
idealized  one-dimensional  model  of  subsidence  which  is  based  on  the 
following  system  of  linear  partial  differential  equations. 


du/dT  =  -ady/dT  (i.l) 

dy/dt  =  u  (->2) 


Here  y(t,T)  is  the  height  above  the  top  surface  of  the  cavity  at  time  t 
of  the  layer  of  material  through  which  the  plane  of  disturbance  passed 
at  time  T;  u(t,T)  is  the  velocity  of  this  layer;  and 


a(t,T)  =  K 


1  +  K(t  -  T 


)]  D 


1  +  t  -  T 


-1 


(1. 


-7  \ 

J) 


is  an  empirical  function,  in  which  the  constant  K  >  1=  By  eliminating 
y(t,T)  between  (i.l)  one  obtains  for  u(t,T)  the  linear  hyperbolic  partial 
differential  equation 

d2u  -  2K  du  (1.4) 

dtdf  "  1  +  K( t-T ) 


Boundary  conditions  for  (1.4)  were  determined  by  the  following  con- 


sidera.ions.  Let  y  be  the  original  density  of  the  undisturbed  overt 


Ui  U.CU  y 


and  7q  the  density  immediately  after  passage  of  the  disturbance,  both 
being  assumed  constant.  Then  conservation  of  matter  at  the  wave  front 


Y  f  1. 1  ^  ImnllpH 
V  w7  ■ - r - - 


7g  dy(T,  T)/dT  =  7q  |_dy(T,  T)/dT  -  u(T,  T)J  (1.5) 

By  (I.l)  -  (l.j)  this  can  be  rewritten  as 

du(T,  T)/dT  -  7S  (K  -  1)  (ye  -  7q)  _1  u(T,  T)  ^.6) 
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On  T  =  0  it  was  assumed,  on  empirical  grounds  that  for  some  constant 
e>0 

u(t,  0)  =  -  €7g  (K  -  1)  (7S  -  70)"1  (1  +  Kt)'2  (1.7) 

Leser  and  Jenike  solved  this  problem  originally  by  numerical  inte¬ 
gration.  However,  as  we  shall  show  in  the  following  sections: 

(i)  Axi  integral  representation  (2.6)  can  be  found  for  u(t,  T)  in 
terms  of  u(t,  t); 

(ii)  u(t,  t)  satisfies  a  Volterra  integral  equation  (2.8)  of  the 
second  kind,  of  convolution  type,  that  depends  on  the  single  parameter 

a  =  (1  -  K"1)  (1  -  7_7  _1)  _1 


o  s 


s.  n 

/  V 


/-I  Q  \ 

\±.v) 


(ill)  lim  u( t,  t)  =  0  or  oo  as  a  =  1  or  a  >  1. 

Xi  o  o 

(iv)  1  =  (l  +  t)  u(t,  t)  =  m  (a)  for  some  m(a),  a  <.1. 

(v)  u(t.  t)  is  a  strictly  increasing  function  of  a  for  fixed  t. 

(vi)  u(t,  t)  is  a  strictly  decreasing  function  of  t  for  fixed  a  such 

that  0  C  a  It  1. 

(vii)  The  wave-front  y(t,  t)  advances  only  a  finite  distance  if  a  <  l. 

The  author  wishes  to  acknowledge  his  indebtedness  to  Dr.  T.  Leser 
for  calling  his  attention  to  this  problem,  ana  to  Dr.  W.  C.  Taylor  for 
suggestions  concerning  its  treatment  by  Laplace  transform  theory. 

REFORMULATION-REDUCTION  TO  AN  INTEGRAL  EQUATION 
If  we  let 

(r„  -  7„) 

(2.1) 


*<*•«  -  -  «SV/;i  »<|.  r  ) 

s 


then  (1.5)  -  (1.7)  become 

2 


>2 
o  w 


8t8T 


8w 


+  x  -  T  OT 


(2.2) 
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8w(T,T) 


=  a  v  (T,  T) 


(2.3) 


where 


w(t,  0)  =  (l  +  t) 


a  -  70  (K  "  i)/K  (rg  -7q)  >  0 


Now  (2.2)  to  (2.4)  imply 


v(t,  T) 


(i  + 1  r 


/o  I.  \ 


(2.5) 


JL 

a  f  ifer)  dr  (a.6) 

/  (1  t  t  -  r)  2 


and  If  we  let 


v(t)  =  w(t,  t) 


then  by  (2.6) 


v(t)  =  — - - o 

(1  +  t  f 


t 

/  v(r)dr 

J0  (i  + 1  -  r; 


\*-  •  i  / 


(2-8) 


Thus  v(t)  is  defined  as  a  solution  of  a  Vol terra  integral  equation  of 
the  second  kind  of  convolution  type.  The  standard  theorem  for  such 
equations  ^ ^  asserts  the  existence  of  v(t),  and  in  fact  yields 


(l  +  t)”  =  v(t)  =  e“°  (1  +  t)' 


Also  by  (2.6) 


(l  +  t)"^  =~  w(t.  T)  v(t,  t)  =  v(t) 


(2.9) 


/  A  1  \ 

(*•10/ 


and,  actually,  w(t,  T)  is  a  strictly  increasing  function  of  T  for  fixed 
t  =  T.  Since  the  positive  function  w(t,  T)  is  dominated  by  v(t),  in 
order  to  provide  a  basis  for  estimating  w  it  will  suffice  to  consider 
only  v(t ) . 

APPLICATION  OF  LAFIACE  TRANSF0RM3 

By  (2.9)  v(t)  is  of  exponential  type.  Accordingly  its  Liprace 
transform 
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_oo 


\r!  „  \ 
»  \°  / 


J 


r  „-st  x 

=  v(x)  ax 


/  —  .V 

13.JJ 


regtu ar  function  of  5  for  Re  s  a.  Also 


f(t)  =  (l+t)‘ 


4  a  V\/m  m/^a^  av\^  *4  4- «  T  o»\1  n»n  -i - -  " - 

4.U  wviauu^>u^  uuu  I/O  ua^xauc  wi  aUDi  Ui  111 


(3.2) 


F(s)  - 


r 

i 

'J  O 


QO 


q  W  4.  \ 

€  J.  \  Is  /  U.  Is 


(3.3) 


is  regular  for  Re  s  >0.  Since 


/  *  0 _ /  \ 


and  since 


-a 


^e  ns;  )'  =  e  /s 


(3.4) 


-St, 

c  r 


-  (e-8F)! 


=  vs  lOr  ne  a  =  oo. 


then  (3.4)  implies,  after  an  integration  by  parts,  that 

~oo 


Ws  ^  =  1  - 


/ 

J. 


-1  -r  , 

’  P  (1r 


which  is  regular  in  the  entire  s  -plane  except  for  a  slit  along  the 


negative  real-axis. 


In  the  region  where  both  V(s)  and  F(s)  exist  the  Laplace  transform 

or  (  ^ .  o  ; 

V(s )  =  F(s)  +  aF(s)  V  (s) 

and  then  (3.5)  and 

V(s)  =  F(s )/  |_1  -  a  F(s)J  (3.6) 

define  the  analytic  continuation  of  V(s)  into  the  slit  s-plane.  To 
determine  the  nature  and  location  of  the  singularities  of  V,  first 
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consider  F.  If  we  move  the  point  s  around  a  closed  path  that  contains 
s  =  0  we  obtain  from  (3.5) 

F  (se2rti)  =  F(s)  +  2jtiseS 
Accordingly 


F(s)  =  H(s)  +  seslog  s  (3-7) 

where  H(s)  is  an  entire  function.  Thus  F(s)  has  a  logarithmic  branch 
point  at  s  =  0.  Since  F  is  real  on  the  positive  real  axis,  H(s)  Is 
real  on  the  entire  real  axis .  Now 


A 


Im  F(s )  =J 


0 


for  s  =  0 


+  ses  for  s^-0 


where  the  upper  (lower)  sign  refers  to  the  upper  (lower)  side  of  the 
slit.  By  the  maximum  principle  for  harmonic  functions  Im  F(s)  Is 
negative  (positive)  on  the  upper  (lower)  half  of  the  s-plane. 
Consequently  the  roots  of 


1  -  aF  (s  )  =  0 

'  o 


(3-8) 


can  lie  only  on  the  positive  real  axis  or  at  s  =0.  Since  F(0)  =  1, 
and  since  for  8=0 


F'(s) 


Oo 


» o+. 


e"“  *  tf  (t)dt  <0 


yo 

(3.8)  has  a  single  simple  root  sq  >  0  if  and  only  if  a  >  1;  and 
8q  =  0  if  and  only  if  a  =  1.  Thus,  if  a  >  1,  then  V(s)  has  a  single 
simple  pole  at  sq  in  addition  to  a  branch  point  at  0;  if  a  <_1  it 
has  only  a  branch  point. 
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The  complex  Inversion  formula  yields 


-  V(t)  =  (2* i)-1 

where  x  >  s  .  An  estimate 
o 

found  by  integrating  (4.5) 
F(s)  =  s"1  -2s 


^-X+i  oo 

j  ets  V(a)  ds 

(3.9) 

-  X- i  oo 

that  can  be  used  to  bound 

J  V(s )  |  can  be 

twice  by  parts  to  obtain 

-2  ,  B 

+  ose 


r 


-4  -r  . 
r  e  dr 


If  the  path  of  integration  of  F(s)  is  taken  to  be  the  half-line 

|  a  j  =  r  =  oo  and  the  circular  arc  r  =  |  a  j  O^L.  |o|=~  0  ^r  «  > 

then 

|F(s)  -  s-1  +  2s“2|  ^6  | s  |  e“2  Is  I  Bln  ^2  +6  | s  J  “2  (3.10) 

By  standard  techniques  and  application  of  (3*10)  the  path  of 
Integration  in  (3*9)  can  be  deformed  Into  a  circle  C  about  sq  (if  a  ^  l) 

and  a  path  that  traverses  both  sides  of  the  negative  real  axiB.  The 

circle  contributes 


(2*1  T1  jc  f  Ta^/I)  =  -  ebotF(so)/aF*(so),  a>  1  (5.11) 

where,  by  (3-5)  and  (3-9) 

F(so)  =  l/a,  aF!(Bo)  =  (l  +  s"1)  (l  -a)  +  a  (3-12) 

The  negative  real  axis  contributes 

f° 

I(t)  =  (Piti)"1  j  e1*  (V“  -  V+  )dx 
'J  _ 


-oo 


where  V~  (V+ )  designates  values  on  the  lower  (upper)  side  of  the  axis 
and  x  =  Re  s  =  s.  Now  write 
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v"  -  V+  =  f“  /(l-aF“)  -  F+/ (l-aF+)  =  (F~  -  F+)/(l-aF~)(l-aF+) 

By  (3.7)  F  =  F  and  F  -  F  =  -  2jtixeA  ,  whence 

r° 


I(t)  =  -  /  xe^1+t'x  |l-aF‘(x)|  ‘2  dx;>  0.  (3.13) 

QO 

If  a  /  1,  then  J  1  -aF  J  /  0  on  00  =  x  ^  0  ,  since  F”(-oo)  ^  0 

and  F~(o)  =  i.  Then  if  we  let 

m(a)  =  min  |l  -  aF- J  ^  0 

t 

we  find 

r  rQ 

0  <l(t)<-  m"2  /  xe^1+t^X  dx  =  m“2(l  ,+  t)-2,  a  f  1.  (3-1^) 

*^CC 

Now  we  can  prove  the  assertions  of  Section  1.,  If  a  >  1,  then 

Sot 

v(t)  -  I(t)  -  e  F(so)/aF»  (sq) 

which  implies  part  of  (iii).  If  a  <.1,  then  v(t)  =  l(t),  and  (3=1^) 

implies  another  part  of  (iii).  In  conjunction  with  v  =  (l  +  t)  ^ 

(3.14)  yields  (iv). 

The  case  a  =  1  remains  to  be  treated.  First  note  that  by  ( 3 -T ) 

F(0)  =  1  implies  H(0)  •>  1.  Now  let 


00 


n 


ms)  =  x  f  /  ^  an  b 
■  1 


where  the  coefficients  a  are  real.  Thei 

n 


oo 


and 


1  -  F-  .  -  JT*n  xn  -  xex  jlog  (-x)  - 
1 


.  1 
JU  • 


J 


r 


li  -  mt  2  .  *2  i 


v”"1  +  ®x  106  (■*) 

i 


2  2x> 
+  it  e  ] 


The  convergence  of  X  at  its  upper  limit  follows  from 


-  f  -  - »  f 

!  .  ,  2  /  x  i  , 


-OO 


_o 


u"^du  =  -  N/log  b  (3.16) 


' -b  x  log~  (-x) 


log  b 


where  0  <.b  <1,  and  N  =  max  (l,  e^  The  convergence  of  I  at 


2  «/  2  2xn  ^ 

-  u(x  e  )  for  x->  -  00. 


Now  we  can  easily  deduce 

llmt=oo  v(t)  =  0  ±£  a  =  X* 


Proof :  Write 


■/“  */ 

y -00  y  -b 


x  e 


(l+t)x  i., 


1  -  F 


-2 


I1  +  I2 


and  suppose  t  >  1.  First,  choose  b  4.1  so  small  that 

o 

I  _  1  y  p  p 

jl  -  F  J  =  0.5x  log  “(-x)  in  I2.  Then  choose  G  ;>  0  and  reduce  b, 
if  necessary,  so  that  -2/log  ’o  4  6.  Then  by  (3.16)  0  c.  Ip  4.  G  for  all 
t  >1.  Next,  if  we  let 


I  x 
'-00 


li 

1 


Ay 


TD 

■L 


then 


_-(l+t  )b 

—  C  ’ 


Pi 


12 


(5-15) 


i 


and  for 


t  >  1  -  b  _1  log  (G/P) 

we  have  0  4  1^  4  G.  Finally ,  for  such  b  and  t,  0  41^  20. 

The  standard  interative  technique  for  solving  (2.8)  would  yield  a 
power  series  in  a  with  coefficients  that  are  positive  functions  of  t . 
Thus  v(t)  is  a  monotone  function  of  a  for  all  t  >  0.  This  and  the 
representation  v(t)  *  l(t)  for  0  4a  =  1,  and  (3.15)  imply  (v)  and 

/..j  \ 

VV-L/. 

Finally,  let  us  consider  the  motion  of  the  wave-front  y  =  y(t,  t). 
By  (1.5),  (2.1),  and  (2.7) 

dy(t.  t)/dt  =  e(K  -l)  7.  7_  (7.  -  7.)  ”2  v(K  t) 


„Kt 


and  if  y(0,  0)  =0,  then 

y(t,  t)  -  e  K"X(K  -  1)  707S(7S  -  70)'2  J  v(s)  ds  (3-17) 
By  (iv),  if  a  /  1. 

t 


-  +  J 


^2 

v(s)  ds  =  m  (a)  , 


which  Implies  (vii). 


^ .  ffc  $US2&L 


J.  H.  GIESE 
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